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Abstract 

In this paper we consider the Grioli-Koiter-Mindlin-Toupin indeterminate couple stress model. The main 
aim is to show that the traction boundary conditions were not yet completely deduced. As it turns out, and 
to our own surprise, restricting the boundary condition framework from the strain gradient models to the 
couple stress model does not reduce to Mindlin’s set of accepted boundary conditions. We present therefore, 
for the first time the complete, consistent set of traction boundary conditions. 
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1 Introduction 

Higher gradient elasticity models are nowadays increasingly used to describe mechanical structures at the micro- 
and nano-scale or to regularize certain ill-posed problems by means of higher gradient contributions [T[ [5]. One 
of the very first among such models is the so called indeterminate couple stress model El II151E] in which 
the higher gradient contributions only enter through gradients on the continuum rotation, i.e. the total elastic 
energy can be written as W(Vu, V(Vrt)) = HA(symVM) + Wc Urv (V(curlw)). 

The question of boundary conditions in higher gradient elasticity models has been a subject of continuous 
attention. The crux of the matter in higher gradient models is the impossibility to vary the test function and 
its gradient independently. A suitable split into tangential and normal parts must always be considered. This 
is well known in general higher gradient models, see e.g. The indeterminate couple stress model has been 

investigated in this respect as well. A first answer has been given by Mindlin and Tiersten |T| as well as Koiter 
[6] who established (correctly) that only 5 geometric and 5 traction boundary conditions can be prescribed due 
to the dependence of the curvature energy only on gradients of rotations. We agree that there are 5 traction 
boundary conditions in the indeterminate couple stress model which may be independently prescribed. However, 
we show in isi uni that the correct traction boundary conditions are not those proposed by Mindlin and Tiersten 
[4] and which are currently used in the literature. Since all papers dealing with the indeterminate couple stress 
model use this incomplete set of boundary conditions we will not refer further to any specific one. 
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2 The indeterminate couple stress model 

We consider a body which occupies a bounded open set SI of the three-dimensional Euclidian space R 3 and 
assume that its boundary dfl is a piecewise smooth surface. An elastic material fills the domain fl C R 3 and 
we refer the motion of the body to rectangular axes Oxi , i = 1,2,3. For vector fields v with components 
Vi £ H 1 (S1), i = 1,2,3, we define Vd = ((V v\) T , (V v 2 ) T , (V^ 3 ) T ) T , while for tensor fields P with the 
rows Pi £ H(div;S2), i = 1,2,3, we define DivP = (divPi,divP 2 ,divP 3 ) T . Equivalently, in index notation: 
(V»)it = Vi,k and (DivP)i = Pij.j ■ In the remainder of the paper, symA and skew A denote the symmetric 
and the skew symmetric part of the matrix A, respectively, tr(A) denotes the trace of the matrix X , ||A|| 
is the Frobenius norm of the matrix X. The identity tensor on R 3x3 will be denoted by 1. We also use the 
operator anti : R 3 —> so(3), so(3) := {A £ R 3x3 \X T = —A}, defined by (anti(u))y = —£ijkVk,^v £ R 3 , 
where Eijk is the totally antisymmetric third order permutation Levi-Civita tensor. We use the curl operator, 
curlu = EijkVkj , V v £ R 3 and denote respectively by • , : and (-, •} a simple and double contraction and the 
scalar product between two tensors of any suitable ordeiQ- Everywhere we adopt the Einstein convention of 
sum over repeated indices if not differently specified. 

The Grioli-Koiter-Mindlin-Toupin isotropic indeterminate couple stress model EHTUSlIS] considers the cur¬ 
vature energy Wc Urv (V(curia)) = ^ || sym Vcurl u\\ 2 + ^ || skewVcurl u || 2 and the classical elastic energy 
kF e (symVw) = /r||symVw|| 2 + j [tr(symVu)] 2 , where /r, A,«i and a 2 are constitutive coefficients. The correct 
and accepted strong form of the Euler-Lagrange equations are 

Div (er — 1 anti(Divm)) + / = 0, equilibrium of forces 

a = P sym vuW / e(sym Vu) = 2 ji symV« + A tr(Vu)l, symmetric Cauchy-stress (2.1) 

fh = Pvcuri«W / curv(VcurlM) = c*i sym(Vcurln) + a 2 skew(Vcurlw), couple stress tensor. 

Note that the couple stress tensor to is a second order and trace free tensor. Having the Euler-Lagrange equation, 
the question of which boundary conditions may be prescribed arises. 


3 The incomplete boundary conditions considered in literature 


We want to stress the fact that in the framework of a complete second gradient theory we can arbitrarily 
prescribe u and the normal derivative of the displacement Vu • n on the Dirichlet boundary T. This means that 
one has 6 independent geometric (or kinematical) boundary conditions that can be assigned on the boundary of 
the considered second gradient medium. Analogously, one can assign 6 traction (or natural) conditions on the 
force (in duality of u) and double force (in duality of Vu • n), respectively, at <9f2 \ T. The situation is slightly 
different in the indeterminate couple stress model since only a certain linear combination of second derivatives, 
i.e. Vcurl it, is controlled. Mindlin and Tiersten [5] concluded that the geometric boundary conditions on 
r C dfl are the five independent conditions 


= u 


(1 — n <3 n) ■ curl u = (1 — n ® n) • curl u° 


(3.2) 


for a given vector function u° at the boundary, where n is the unit normal vector on dfl and (8> denotes the 
dyadic product of two vectors. The latter condition, in fact, prescribes only the tangential component of curlu. 
Therefore, one may prescribe only 5 independent boundary conditions. 

The possible traction boundary conditions on the remaining boundary dfl \ F given first by Mindlin and 
Tiersten [3] are 


{ (cr — A anti(DivTO-)) • n — x V[(n, (sym?n) • n)]} 

(1 — n <g> n) ■ fh ■ n 

for prescribed vector functions t and g at the boundary, where (•, •) denotes the scalar product of two vectors. 
Mindlin and Tiersten [3] have correctly concluded that the maximal number of independent traction boundary 
conditions is also 5. The same conclusion has been arrived at by Koiter |6|. These traction boundary conditions 
(El have been rederived again and again. However, they are erroneous. 

For example, (T.v)j — AijVj , A . B — A^j Bji , (C ■B'jijfc — CijpBpfc, (O . B'ji — Pijp^pji — V'W — 

Vi U'i , (.'1 ■ B j — AijBij etc. 


_ =t, 
dn\r 

„ = (1 — n ® n) • g, 

an\r 


(3.3) 
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4 The correct boundary conditions in the indeterminate couple stress 
model 


The prescribed traction boundary conditions (13.31) proposed by Mindlin and Tiersten |3j do not remain indepen¬ 
dent, in the sense that g leads to a further energetic conjugate, besides t, of u. From this reason and looking back 
to the clear and correct boundary conditions considered in the more general second gradient elasticity model, in 
order to prescribe independent geometric boundary conditions and their corresponding completely independent 
energetic conjugate (traction boundary conditions), we have to prescribe u and (1 — n 0n) • (Vu • n). Let us 
remark that prescribing u| = u° and (1 — n ® ro) • (Vw • n) | r = (1 — n 0 n) • Vu° • n is fully equivalent with 
prescribing it| r = u° and (1 — n® n) • curlu| r = (1 — n® n) • curlu 0 , which is (13.21) . However, in the formulation 
of the principle of virtual power, the energetic conjugate of (1 — n ® n) ■ curl it is not equal to the energetic 
conjugate of (1 — n 0 n) • (Vu • n). 

Using the principle of virtual power proposed by Mindlin and Tiersten [4j Eq. (5.13)], but now suitably 
applying the surface divergence theorem dunum we arrive at the following traction boundary conditions on 

dfl \ T 


| ^er —1 anti(Div?n)^ • n — x V[(n, (synun) • n)] 

— iV [ anti ( (1 — n ® n) ■ fh ■ n ) ■ (1 — n 0 n) ] : (1 — n ® n) j 


an\r 


= t, 


(4.4) 


(1 — n 0 n) • anti[(l — n 0 n) ■ m ■ n] ■ n _ = (1 — n®n)-o, 

an\r 


together with the traction boundary conditions on <9T 

{([anti[(l — n 0 n) ■ fh ■ n]] + — [anti[(l — n 0 n) ■ fh ■ n]] - ) • u} 


a r 


(4.5) 


where t and g are prescribed vector functions on \ T, while n is a prescribed vector function on <9T. Here, 
v is a vector tangential to the surface T and which is orthogonal to its boundary <9T. The term [anti[(l — n 0 
n) ■ m • n]] + — [anti[(l — n 0 n)-m ■ n]] - measures the discontinuity of anti[(l — n ® n) ■ fh ■ n] across <9T. 

Comparing (13.31) and (14.41) . we remark that in the Mindlin and Tiersten formulation (13.3IL it remains a 
missing boundary term — -i V [ anti ((1— n® n) • fh • n) • (1 — n®n)] : (1 — n 0 n) which also performs work 
against u. On the other hand, we show in [9] that when the higher gradient contributions only enter through 
gradients on the continuum rotation, i.e., Vcurlu, the independent traction boundary conditions which are 
coming from the representation in terms of a particular case of second gradient elasticity model written with 
third order moment tensors coincide with our novel traction boundary conditions (14.41) and (14.51) , and not with 
the traction boundary conditions (13.31) proposed by Mindlin and Tiersten. 

Our renewed interest in traction boundary conditions in the indeterminate couple stress model was triggered 
by the controversial papers mm ■ There, the authors have made far reaching claims on the possible anti¬ 
symmetric nature of the second order couple stress tensor fh. Their reasoning is based on physically plausible 
assumptions similar to a Cosserat or micromorphic theory m which led them to require a total split of the 
effect of force and moment tensors in (13.31) . In (13.31) . this can be achieved if and only if fh is skew-symmetric and 
this constitutes the essence of their claim. However, since (El is incomplete, their conclusion is misleading, see 
also m ■ The couple stress tensor in the indeterminate couple stress theory is not necessarily skew-symmetric! 
Quite to the contrary, the couple stress tensor may be chosen to be symmetric mm- 
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